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3 CHE  COMMENTS  AND  FROSFECTITOS  ON  DTE  DIFFUBICN  EXPERIMENTS 
(Unpublished  Technical  Note) 

by  T&kashl  I  oh  lye 


1.  INTRODUCTION 

This  note  is  c  result  of  the  author's  efforts  In  groping  for 
sons  interpretation  of  dye  diffusion  data  accumulated  at  the  Lamont 
Geological  Observatory  since  1961*  It  is  also  the  purpose  of  this 
note  to  give  sene  ideas  on  desirable  future  experiments  and  to  sti¬ 
mulate  theoretical  study  on  diffusion  in  the  ocean.  Most  of  the  data 
referred  to  are  described  in  the  unpublished  Technical  Paper  "Dye 
Diffusion  Experiments  in  the  New  York  Bight”  by  Cost  in,  Davis,  Gerard 
and  Kata  (1963),  cited  as  ENT.  .  . 

2.  (RGBS  FEATURES  OF  DTE  PATCHES 

As  described  in  ENT  and  by  Iohiye  (1962a),  visual,  photographic 
and  fluorometrlo  patterns  of  dye  patches  show  a  tendency  of  elongation 
in  the  direction  of  surfaoe  currents  (or  winds  when  they  are  suffici¬ 
ently  strong).  There  are  two  possible  causes  for  such  elongation. 

One  cause  is  that  the  turbulanoe  energy  (defined  as  the  mean  square 
fluctuation  velocities)  is  more  Intense  in  the  direction  of  the  mean 
current  than  acroea  the  mean  current.  The  other  is  that  the  mean 
current  has  a  vertical  shear,  usually  with  a  maximal  velocity  at  the 
surfaoe,  and  thus  the  dye  concentration  at  each  level  spreads  with 
different  speeds  in  the  direction  of  the  mean  current.  In  an  equili- 
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brlvm  state,  it  is  believed  that  turbulence  tends  to  be  isotropic. 

Hew  ever,  there  is  no  reason  to  generalise  an  assumption  of  isotro  pi- 
city  to  the  surface  current  generated  by  transitional  wind  stresses. 
Also  there  is  experimental  evidence  that  large  eddies  in  a  shear  flow 
are  elongated  in  the  direction  of  the  flow  (Townsend,  1956).  There¬ 
fore,  the  first  cause  might  be  a  possibility,  but  experimental  veri¬ 
fication  in  the  sea  and  theoretical  models  of  such  anisotropic  tur¬ 
bulence  are  yet  to  come. 

The  second  cause  can  be  analysed  by  considering  the  diffusion 
equation  of  the  type* 


in  which  S  is  the  concentration  of  dye,  0  is  the  mean  horizontal 
velocity  varying  in  the  vertical  direction  Z,  and  Ag,  Ay  and  A^ 
are  diffuaivity  in  x,  y  and  s  directions,  respectively.  Application 
of  dye  technique  to  determine  flow  rate  and  diffusion  in  a  pipe  or 
an  open  channel  was  already  discussed  by  Taylor  and  Elder,  using  the 
diffusion  equation  (1).  See  the  author's  review  (lohiye,  1962b) 
pp.  U2~2l6.  However,  the  flow  in  a  pipe  or  channel  la  much  simpler 
than  in  the  sea.  There  is  no  reason  to  believe  that  the  results  of 
experiments  In  a  small  hydraulio  fltns  a re  ocatplstsly  valid  In  the 
sea. 

Another  interesting  problem  concerning  the  mean  currant  is  the 
curvature  of  dye  patterns  discussed  by  Ooetln,  Davis,  Qerard  and  Eats 
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In  "Ogre  Diffusion  Experiments  In  the  New  York  Bight"  (1963).  They 
described  many  examples  of  dye  patches  which  showed  counter-clockwise 
movement  of  surface  water  relatively  to  the  deep  water*  They  estimated 
this  movement  from  observations  indicating  that  the  maxlsnaa  concentra¬ 
tion  contained  by  the  surface  water  moved  in  most  eases  cyclonically 
against  the  tail  part  in  a  curved  pattern.  They  explained  such  curved 
patterns  from  the  Ekman  spiral  in  the  upper  layer  of  the  sea,  since 
effects  of  change  in  tidal  currents  and  wind  directions  with  time  seem 
to  be  too  irregular  In  all  series  of  experiments  to  produce  systematic 
anticyclonic  (from  the  surface  to  the  bottom)  curvature  of  the  ho  do- 
graph.  Their  discovery  is  quite  interesting  In  testing  Ekman*  s  theory 
of  wind  drift,  although  possibilities  of  effects  of  tidal  currents  and 
wind  shifts  on  the  curved  patterns  are  still  not  quits  eliminated. 

There  are  two  approaches  to  studying  oceanic  turbulence,  using 
curved  patterns  of  dye  patches,  whan  vertical  distributions  of  currents 
and  dye  concentration  are  known*  One  approach  is  to  compare  observed 
patterns  of  dye  pat  dies  with  a  solution  of  the  diffusion  equation  and 
determine)  pavadietera  such  as  eddy  diffuaivity.  One  such  equation  is 
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In  which  vertical  diffusion  is  neglected,  but  x  and  y,  components  of 
me ?in  current,  are  included,  taking  Into  account  spiral  structures  of 
the  wind  drift*  A  solution  of  equation  (2)  for  the  initial  source, 
which  is  expressed  by  v\'x,  j,  2)  Is  given  by* 
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When  the  vertical  mixing  is  considered)  the  sofution  of  a  diffusion 
equation  becomes  more  complicated  than  equation  (3),  but  it  can  be 
obtained  either  numerically  or  by  successive  approximation  from  (3). 
However,  there  is  no  need  to  determine  the  exact  solution,  because 
the  diffusion  equation  itself  is  an  approximation.  Further,  the 
quantities  appearing  in  the  diffusion  equation  can  be  measured  in  the 
sea  only  with  limited  accuracy.  In  most  cases  of  the  experiments, 
only  the  solution  for  a  vertical  line  source,  which  is  given  by  F(z), 
is  sufficient  for  determining  order  of  magnitudes  of  eddy  dif  fusivity. 
The  concentration  in  aerial  photographs  of  dye  patches  is  represented 
by  an  integral: 


V*  S1^’  3  ,  £,<r)  e  ’■kzdz  m 

C 

in  which  k  is  a  vertical  extinction  coefficient  for  visibility  of  dye 
in  the  water,  and  »-ax±s  is  taken  downwards  from  the  surface.  The 
coefficient  k  is  to  be  determined  from  experiments,  in  which  dye 
patches  of  known  concentration  at  different  depths  must  be  photo¬ 
graphed  from  an  airplane. 

Another  approach  to  oceanic  turbulence  by  using  curved  patterns 
is  to  determine  vertical  eddy  viscosity  in  the  upper  layer.  Suppose 
that  dye  is  initially  distributed  uniformly  along  a  vertical  axis. 


-*• 


If  there  is  a  horizontal  current  whose  direction  and  speed  are  vari¬ 
able  with  depth,  the  initial  dyed  filament  of  water  is  distorted. 

If  there  is  no  diffusion,  projection  of  the  filament  on  a  horizontal 
plane  represents  a  hodograph  (of  the  trajectory),  and  it  forms  a 
spiral  for  wind  drifts  corresponding  to  Ekman's  theory.  The  projected 
curve  can  be  expressed  by  a  position  vector  j\_  ”  (  ^  (5) 

defined  byi 


in  which  the  initial  values  of  4?  and  ^  are  taken  at  the  origin  of 
coordinates.  The  curvature  fl~l of  the  projected  curve  of  the  fila¬ 
ment  (called  hereafter  simply  the  dyed  curve)  is  given  by  a  function 
of  the  distance  s  measured  from  an  arbitrary  point  along  the  curve, 
such  asi 
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according  to  a  theorem  of  differential  geometry,  where « 


d-6^  -  Jg*  + 


(8) 


In  order  to  express  the  curvature  ty  functions  of  a  parameter  Z 
Instead  of  s,  the  following  relations  are  used* 


d2^  fclz\\  d d*£ 
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in  which  the  subscript  z  indicates  differentiation  with  z.  Substi¬ 
tution  of  (9b)  and  (9c)  into  (9a)  yields i 

r  (X)-1  =  1 7zh*r  hfa  i  ( s£  *  ?*•) ■  v*  (»> 

In  a  stationary  and  horizontally  uniform  current,  £?  and  A? 
becomes t 
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and  equation  (10)  can  be  expressed  in  terms  of  horizontal  velocity. 
When  the  current  is  generated  from  balance  of  Coriolis'  force  and 
shearing  stresses,  as  in  Elman*  s  wind  drift,  the  equations  of  motion 


are  i 
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in  which  f  is  the  Coriolis*  ooeffiolent  and  k(n)  Is  the  eddy  viscosity 
variable  with  depth.  Substituting  relations  of  j^^and  l/^ot  (12 
a  !t  b)  into  (10),  we  have! 


K(tO  -  -  . 


(33) 
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Equation  (13)  Indicates  that  the  eddy  viscosity  at  each  level  can 
be  determined  when  the  curvature  of  the  dyed  curve  and  velocity  and 
shear  of  the  mean  currents  are  known*  It  is  noted  that  the  eddy 
viscosity  is  proportional  to  the  radius  of  curvature*  When  the 
longer  central  line  of  an  elongated  pattern  is  interpreted  as  the 
hodograph  of  the  wind-driven  current,  the  example  of  the  patch  in 
figure  9  of  EWT  gives  the  following  numerical  valuest 

4  -  \Q~^  Csec'O  ;  R  ^  I O5"  Ccw) 

U  /v.  V*  ^  10  (c  **  /sec)  •  Ujr  |  o'1  rsec'1} 

t  --  I  o  ^  (Sec) 

Equation  (13)  gives  the  averaged  value  with  depth  of  K  equaling 
100(cm^/8ec),  which  is  not  unreasonable  compared  with  other  estimates 
of  vertical  eddy  viscosity*  In  more  general  cases  of  currents  which 
are  nan-stationary  and  driven  by  pressure  gradients  as  well  as  winds, 
an  equation  for  k(z)  similar  to  (13)  can  be  obtained  fay  substituting 
equations  of  motion  into  (10),  beoause  the  tern  containing  k*  (») 
vanishes  in  the  right-hand  side  of  (10).  However,  since  the  resulting 
relation  has  so  many  terms  which  are  not  easily  estimated  from  obser¬ 
vation,  such  elaboration  is  hardly  justifiable. 

3*  APPLICATION  CF  ISOTROPIC  DIFFUSION  MODELS 

Data  on  dye  diffusion  were  very  poor  in  their  quality  a  few  years 
ago,  since  there  was  no  continuously  recording  type  fluoroneter  and 


even  good  serial  photographs  were  scarce.  Hi  order  to  obtain  a 
structure  of  turbulence  from  such  data,  the  author  (Ichlye,  1959) 
determined  change  of  areas  with  $re  concentration  above  a  threshold 
value  using  visual  observations  of  dye  patches  and  applied  Taylor* s 
theory  of  diffusion  by  discontinuous  movements  and  a  semi- empirical 
theory  of  Joseph  and  Sandner,  assuming  that  the  diffusion  of  dye 
patches  is  horizontally  isotropio.  Okubo  (196?)  elaborated  this 
approach  and  tried  to  verify  various  diffusion  models  (isotropio) 
by  curve-fitting  techniques.  However,  such  a  method  of  verifica¬ 
tion  is  almost  meaningless,  since  any  curve  may  fit  (or  deviate  from) 
few,  scattered  points  of  observed  data  of  dubious  quality.  (This  doss 
not  mean  that  theoretical  study  on  diffusion  process  is  meaningless, 
though. ) 

Mary  aerial  photographs  of  high  quality,  combined  with  concen¬ 
tration  contours  determined  with  fluoro meters,  indicate  that  details 
of  dye  patches  with  elongated,  finger-like  or  striated  shapes  are 
quite  different  from  the  features  predicted  by  a  simple  theory  of 
isotropio  diffusion.  Also,  as  indicated  In  ENT,  other  environmental 
factors  suggest  that  the  shape  of  these  patches  is  dependent  on  ad- 
v active  currents  and  vertioal  diffusion,  as  discussed  in  the  previous 
section,  or  on  instability  of  shearing  flow,  as  will  be  discussed  in 
the  next  section.  It  is  the  zeroth  order  approximation  to  the  diffu¬ 
sion  process  to  apply  results  of  conventional  theories  mi  turbulent 
diffusion  to  dye  patterns  assumed  as  oiroumsynaetrioal  and  derive 
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turbulence  characteristics  by  curve  fitting.  The  next  approximation 
must  be  such  that  the  patterns  of  dye  patches  may  be  correlated  with 
direct  measurements  of  mean  currents,  their  fluctuations,  stability 
and  other  environmental  data. 

It.  STRIATED  STRUCTURES  CF  DYE  PATCHES 

Mary  examples  of  dye  patches  described  in  ENY  showed  striated 
patterns,  in  which  a  rather  regular  array  of  furrows  is,  in  most  cases, 
almost  parallel  to  the  wind  direction.  In  some  cases,  when  the  dye 
patch  was  elongated  and  curved,  as  described  in  section  2,  the  outlook 
of  the  striated  pattern  superimposed  on  the  patch  is  similar  to  the 
pattern  of  flows  caused  by  roll-type  instability  in  rotating  tank 
experiments  of  Faller  (1963).  Howsver,  such  similarity  might  be  only 
superficial,  because  scales  of  streaks  and  time  period  of  their  deve¬ 
lopment  are  eo  mall  In  the  natural  water  that  effects  of  Coriolis' 
terms  seem  to  be  unimportant  to  generation  of  such  streaks.  But 
Coriolis'  force  might  be  effective  when  the  wave  lengths  (or  interval 
between  the  furrows)  of  streaks  are  large  and  the  striated  pattern  is 
maintained  for  a  long  period.  To  check  such  effects  theoretically  is 
s  challenging  problem  for  geophysical  hydrodynami cists. 

Mechanism  of  generation  of  striated  patterns  seems  to  be  quite 
intriguing.  If  it  is  assumed  that  the  furrows  are  due  to  a  roll-type 
oellular  motion  caused  by  instability  in  a  uniform  surface  current 
with  vertical  shear,  such  motion  must  be  generated  by  non- linear  inter¬ 
action  between  the  mean  current  and  perturbations.  This  can  be  easily 


seen  from  linearised  perturbation  equations  for  a  shear  flow  discussed 
by  Lin  If  the  wave  velocity  of  perturbations  have  a  component 

moving  parallel  to  the  mean  current,  stability  of  the  current  can  be 
discussed  as  a  linear  perturbation  problem  as  treated  by  Lin.  However, 
if  the  perturbations  are  variable  only  in  the  directions  perpendicular 
to  the  mean  current,  the  linearised  perturbation  equations  yield  only 
gravity  waves  (with  damping)  with  crests  parallel  to  the  mean  current* 
Therefore,  in  order  to  determine  the  stability  of  the  shear  current 
for  the  long  perturbations  parallel  to  the  current  as  in  cellular 
motion  which  is  assumed  to  cause  furrows,  non-linearity  must  be  con¬ 
sidered*  This  is  a  reason  that  mathematical  formulation  for  flows 
responsible  for  Lagmuir  streaks  is  not  achieved  yet. 

Welander  (1963)  proposed  a  speculation  on  mechanism  of  natural 
streaks  on  the  sea  surface.  He  assumed  that  material  forming  a  sur¬ 
face  film  of  the  sea  might  be  concentrated  by  chance  along  a  band  in 
the  direction  of  the  wind,  even  when  the  wind  is  uniform.  Such  material 
causes  damping  of  capillary  waves  in  the  band.  The  wind  becomes 
stronger  along  the  band  due  to  decrease  of  surface  friction  than  on 
both  sides  of  it*  Such  difference  in  friction  will  also  create  a 
secondary  cellular  circulation  perpendicular  to  the  mean  wind,  and 
the  air  close  to  the  surfaoe  flows  toward  the  band*  Such  air  motion 
in  turn  will  produce  cellular  motion  of  the  water* 

Welander1 s  hypothesis  seems  to  be  novel,  because  the  mechanism  is 
sought  in  the  air  motion  rather  than  in  the  water  motion,  and  also  the 


Importance  of  a  surface  film  is  assumed*  However,  his  hypothesis 
cannot  explain  the  regularity  of  the  intervals  of  furrows  observed  in 
dye  patches.  In  order  to  test  his  speculation  and,  at  the  same  time, 
to  determine  fine  structures  of  the  streaks,  three  measurements,  which 
are  feasible  with  relative  ease,  can  be  suggested  in  connection  with 
observations  of  dye  patches.  The  first  is  to  measure  fine  structures 
of  wind  over  the  sea  by  using  a  smoke  signal.  The  seeond  is  to  take 
underwater  pictures  (movie  or  still)  of  rear  views  of  striated  patches, 
either  by  a  sklndiver  or  with  an  underwater  TV  camera.  The  third  is 
to  measure  in  situ  surface  tension  of  the  water  or  to  change  it  arti¬ 
ficially  by  adding  chemicals  to  the  sea  surface. 
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